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ABSTRACT 

We present the detection of a radio-emitting jet from the black-hole candidate and 
X-ray binary source Cygnus X-l. Evidence of a bright core with a slightly extended 
structure was found on milliarcsecond resolution observations with the VLBA at 15.4 
GHz. Later observations with the VLBA (and including the phased up VLA) at 8.4 
GHz show an extended jet-like feature extending to ~ 15 mas from a core region, with 
an opening angle of < 2°. In addition, lower resolution MERLIN observations at 5 
GHz show that the source has < 10 per cent linear polarization. The source was in the 
low/hard X-ray state during the observations, and the results confirm the existence of 
persistent radio emission from an unresolved core and a variable relativistic (> 0.6 c) 
jet during this state. 
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^ ■ 1 INTRODUCTION 



3 



Of the ~250 currently known X-ray binaries, approximately 
50 have detectable radio emission and a dozen or so of these 
have been found to have radio jets. Extreme examples of 
the jet sources are GRS 1915+105 (Mirabel & Rodriguez 
1994; Fender et al. 1999) and GRO J 1655-40 (Hjellming & 
Rupen 1995; Tingay et al. 1995), both exhibiting apparent 
superluminal motion. Such jet sources often show strong ra- 
dio flux variability (see Fender, Burnell & Waltman 1997 
for a review). However, some sources show relatively little 
variation, such as Cygnus X-l and LS 5039. The persistently 
radio-emitting (20-40 mjy) LS 5039, which is at a similar 
distance, has been shown to have a milliarcsecond-scale jet 
(Parades et al. 2000). 

Cygnus X-l (HDE 226868, V1357 Cygni) comprises a 
supergiant secondary (spectral type 09.7 lab) with mass be- 
tween 20-33 Mq , together with a compact primary which is 
a strong black- hole candidate with mass between 7-16 Mq. 
The orbital period has been detected in optical, UV, X-ray 
and radio wavelengths (Pooley, Fender & Brocksopp 1999; 
Brocksopp et al. 1999b). The orbital period is <~5.6 days or- 
bital radius of 0.2 AU is derived; at an assumed distance of 
2 kpc (Gierlihski et al. 1999) this corresponds to an angular 



scale of 0.1 milliarcsecond (mas). An orbital ephemeris has 
been derived from optical radial velocity data (LaSala et al. 
1998) and compared with other wavelength data (Brocksopp 
et al. 1999b). 

The source has been fairly stable in the low/hard X- 
ray state, at around 15 mjy at cm wavelengths and with 
an extremely flat radio spectrum. The radio emission was 
first observed following a 'turn on' corresponding to a X- 
ray state change in 1972. The 'turn on' was probably the 
recovery from a high/soft state during which time the jet 
had been quenched as Zhang et al. (1997) show that the 
radio emission was similarly suppressed during the high/soft 
X-ray state in 1996. This is behaviour similar to that shown 
by GX 339-4 during 1998 (Fender et al. 1999b). 

At 15 GHz, a quasi-sinusoidal 3 mjy semi-amplitude 
modulation is observed superimposed on the 14 mjy mean 
level. The minimum occurs at superior conjunction of the 
compact body, implying that the radio-emitting region is 
associated with this compact object (Pooley et al. 1999). 
More recently a further modulation period of 142 days has 
been detected (Brocksopp et al. 1999b), explained as pre- 
cession or warping of the disk. The source displays a flat 
spectrum between 2 and 220 GHz, with an average spectral 
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Figure 1. The naturally weighted VLBA image of Cygnus X-l 
from 1997 March 28, at 15.4 GHz has a lowest contour of 0.4 mjy 
beam -1 and a peak flux density of 5.0 mjy beam" 1 . As with all 
images in this paper the contouring represents factors of two. 



index \a\ < 0.15 (Fender et al. 2000). A small radio flare co- 
incident with an X-ray flare in 1975 further demonstrated an 
X-ray/radio correlation (Hjellming & Han 1995). Both the 
observed flat spectrum and the orbital modulation fit into 
the slowed adiabatically expanding conical radio jet model 
(Hjellming & Johnston 1988). It is not clear how common 
the small flares (less than ~30 mjy) are or exactly how the 
X-ray state changes affect the radio emission (Brocksopp et 
al. 1999b). 

An accretion disk is essential in models of the X-ray 
variability. The source has high/soft, low/hard and interme- 
diate spectral states, spending most time in the low/hard 
X-ray state. The X-ray state changes can be modelled us- 
ing a variable accretion rate; when the accretion rate is low 
an advection dominated accretion flow (ADAF) is postu- 
lated (Esin et al. 1998), giving a low radiative efficiency, as 
expected for such underluminous quiescent sources. The ex- 
istence of a jet perpendicular to the disk implies that the 
outflow process may be important in terms of mass and an- 
gular momentum loss from the disk. 

Previous radio maps on both Very Large Array (VLA) 
(Marti et al. 1996) and MERLIN (Newell 1997) size scales (5 
arcsecond and 50 mas beam FWHM respectively) showed no 
conclusive evidence for resolved or extended emission. The 
VLA image showed emission regions to the north and south 
which could be lobes associated with the source at an extent 
of ~5 arcminutes. 



2 VLBI AND MERLIN OBSERVATIONS 

For over a decade Cygnus X-l has been employed as an 
astrometric reference source (Lestrade et al. 1995). Varia- 
tions in the visibility amplitude were observed on the very 



longest transatlantic baselines (Lestrade, priv. comm.) sug- 
gesting resolved source structure on the milliarcsecond scale. 
To investigate the milliarcsecond source structure in more 
detail, we obtained full track, 15.4 GHz Very Long Baseline 
Array (VLBA) observations of the source in 1997. These 
were subsequently followed up in 1998 with multi-epoch and 
multi-frequency 8.4 and 15.4 GHz observations. The initial 
15.4 GHz observations and the later 8.4 GHz multi-epoch 
observations are discussed in this paper. The multi-epoch 
15.4 GHz observations will be presented in a separate paper 
(de La Force et al. in prep) . In this paper we also present the 
results of an analysis of archive MERLIN data, investigating 
the linear polarization properties of Cygnus X-l. 

2.1 Initial 15.4 GHz VLBA observations 

On 1997 March 28 the VLBA interferometer network was 
used to observe a single epoch of Cygnus X-l at 15.4 GHz 
while the source was in the low/hard X-ray state. The run 
was made at 10:15 - 19:45 UT (centred on Julian date 
2450901.06), recording at a sustained rate of 128 Mbits per 
second. Only one hand of polarization was recorded in Left- 
hand Circular Polarization (LCP) providing a total band- 
width of 64 MHz. The observations were made in phase- 
reference mode with a typical cycle time of 2.0 minutes on 
the phase-calibrator and 2.0 minutes on the target, Cygnus 
X-l. The phase-calibrator used throughout these observa- 
tions was B1955+335, a bright (Sisghz ~300 mjy), com- 
pact radio source (with an accurately measured astromet- 
ric ICRF position) lying within 1.1 degrees of Cygnus X-l. 
The data were correlated at the NRAO correlator in So- 
corro, NM, USA. The coordinates chosen to correlate the 
data were those given by the most recent published radio 
position, 19 h 58 m 21:680 +35°12'05'.'887 (epoch 1986, J2000) 
Lestrade et al. (1995) . With a total of 8 hours observing time 
excellent ttv-coverage of the target source was obtained. 

The initial analysis of the data was performed with the 
NRAO MPS package (e.g. Fomalont 1981). The visibility 
amplitudes were calibrated using the system temperatures 
and gain information provided by each of the telescopes. 
Occasional observations of 3C345 were used to align the in- 
dividual 8 MHz sub-bands (IFs) that spanned the total 64 
MHz bandwidth. The antenna residual delays, fringe rates 
and gain corrections were determined from the phase cali- 
brator and interpolated and applied to the target, Cygnus 
X-l. 

Unfortunately, the positions we supplied to the correla- 
tor did not take into account the proper motion of Cygnus 
X-l. Thus our first map made from the calibrated Cygnus X- 
1 data showed no radio emission in the centre of the map but 
obvious side-lobe structure from a source apparently located 
outside the field was clearly visible. A wide field image of the 
region (for which no averaging of the data were performed) 
revealed this source to be Cygnus X-l, offset 46 mas to the 
east and 84 mas to the south of the nominal phase-centre. 
Our position, after correcting for the measured proper mo- 
tion of the source (provided by Lestrade, priv. comm.) is 
within 1.5 mas of the original astrometric position (Lestrade 
et al. 1995). The refined proper motion is 8.7 ± 0.2 mas yr _1 
in position angle (PA) -151 ± 1°. 

The Cygnus X-l data were Fourier transformed and 
then CLEANed using the AIPS task IMAGR which allows 
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variable data weighting using the robust parameter (Briggs 
1995). No self-calibration of the data was performed. A nat- 
urally weighted (robust 5) image is presented in Figure 1. 
The elliptical Gaussian restoring beam measures 1.13 x 0.63 
mas in PA 12.9°. 

The source is clearly detected, and is plotted relative 
to the epoch 1986 coordinates chosen for correlation. The 
structure is confused and ambiguous, this is likely due to 
errors in the phase referencing process. Unfortunately the 
source was too weak for reliable self-calibration with this 
observing array. An attempt to model fit these data as a sin- 
gle component with DIFMAP (Pearson et al. 1994) suggested 
an extension in PA —15° ± 5 (Stirling, Spencer & Garrett 
1998). It seems likely that the radio emission represents an 
unresolved radio-jet. 



2.2 Follow up VLBA multi-epoch 8.4 GHz 
observations 

2.2.1 Observations and data reduction 

Further high resolution observations of Cygnus X-l were 
made in August 1998 with the VLBA and VLA. Three 
epochs of observations were obtained every 2 days to cover a 
large fraction of the period (5.6 days) of the binary system. 
The runs were made at 04:19 - 08:55 UT on the 10th, (epoch 
A) at 04:20 - 08:30 UT on the 12th (epoch B) and at 04:20 
- 08:30 UT on the 14th (epoch C). The mid point of each 
observation was on Julian date 2451035.76 (A), 2451037.76 
(B) and 2451039.76 (C). Observations were made at 15.4 
and 8.4 GHz, the frequencies being alternated each hour. 
Only the 8.4 GHz data are considered in this paper. 

Since the source is relatively weak, the phased VLA was 
also included as one element in the array. A data rate of 256 
Mbits per second was used (recording LCP in 8 x 8 MHz 
sub-bands, 2-bit sampled at the Nyquist rate). The Cygnus 
X-l data were correlated using a position of 19 h 58 m 21! 67634 
+35°12'05'.'7940, J2000 (derived from our earlier 15.4 GHz 
observations but with the source proper motion now in- 
cluded). The data analysis followed the standard procedure 
described earlier in Section 2.1. 



2.2.2 Images and component fluxes 

Due to thunder storms and bad weather at the VLA the 
third epoch had generally noisier data on the most sensitive 
VLA-VLBA baselines. The three naturally weighted images 
at 8.4 GHz have been self-calibrated in MPS (using the ex- 
tremely sensitive 'phased' VLA as a reference antenna) and 
are presented in Figure ti with identical contouring. Here the 
radio emission from a jet to the north-west is much clearer, 
extending out to ~15 mas from the core. The unresolved 
core likely represents the base of the visible jet and its flux 
density variations with orbital phase are given in Table fU 

There is some correlation between orbital ephemeris 
(from Brocksopp et al. 1999a) and the flux densities mea- 
sured from our maps, as the unresolved core flux density has 
a maximum around inferior conjunction (see Table uh and 
the orbital phase modulation of the radio flux does indeed 
imply the radio flux variation should originate in or near the 
core, as observed. The flare preceding our observations by 
around 15 days as shown in Brocksopp et al. (1999b) may 



complicate the integrated flux density and period relation- 
ship found by Pooley et al. (1999). 



2.3 MERLIN archival data 

As all of the VLBA observations were made using a single 
hand of polarization, it was impossible to examine the linear 
polarization emitted from the jet. As MERLIN observes in 
all four polarization correlations as standard, linearly polar- 
ized emission can be imaged. Only the total intensity (Stokes 
I) image was presented by Newell (1997). 

These MERLIN data, observed on 1995 July 19 were 
retrieved from the archives and calibrated in full polariza- 
tion using the MERLIN pipeline (Thomasson et al. 1993). 
The 5 GHz observations were made in the standard MER- 
LIN continuum setup over a bandwidth of 16 MHz. Also 
observed were a phase reference source, 1951+355; a flux 
and polarization PA calibrator, 3C286; and a point source 
calibrator, OQ208. The total intensity map of Cygnus X-l 
was convolved with a Gaussian beam of 58 x 37 mas in PA 
-81.2° and found to be a point source with a peak flux den- 
sity of 11.9 ± 0.6 mjy beam -1 in agreement with Newell 
(1997). No linear polarization was detected associated with 
the source to an upper limit of 0.9 mjy beam -1 , ruling out 
a simple, very ordered field structure. 



3 ANALYSIS OF THESE AUGUST 1998 DATA 

3.1 Emission mechanism 

We can place a lower limit on the brightness temperature of 
the unresolved radio emission simply by using the Rayleigh- 
Jeans approximation. Using a flux density of 6 mjy and a 
solid angle equivalent to the uniform beam (epoch B) we de- 
rive a minimum temperature of 10 r K at 8.4 GHz, inferring 
that the radiation is non-thermal, and is not for example due 
to free-free emission from gas at 10 4 K. The observed spec- 
tral index and the detection of polarization in a number of 
jet sources suggests that the radiation is via the synchrotron 
mechanism. A spectral break is not observed in the radio or 
mm regimes, implying that the emitting electrons do not 
suffer significant radiative losses (Fender et al. 2000). The 
flat spectral index supports the non-thermal interpretation, 
and is often attributed to multiple expanding synchrotron- 
emitting regions within the beam when observed in the jets 
from quasars. 

3.2 Doppler boosting 

Following the method of Mirabel & Rodriguez (1994), under 
the assumption of an intrinsically symmetric ejection, then 
for bulk motions at a jet velocity the observed ratio of 
the approaching jet components flux density, Sa, PP , to the 
receding jet components flux density, SVcc, (assuming op- 
tically thin components or non-overlapping optically thick 
components) is given by 



Japp 



1 + cos 



Srcc {l-0COSi) (1) 

with S„ oc v +OL ; k = 2 for a continuous jet and k = 3 for 
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Figure 2. VLBA and phased VLA images of Cygnus X-l from August, 1998, at 8 GHz; lowest contour 0.1 mjy. Epoch A is on the left 
with peak flux density 8.7 mjy beam -1 and convolved with a Gaussian beam 2.88 X 1.19 mas in PA -10.9°, epoch B in the middle with 
peak flux density 9.0 mjy beam - and convolved with a Gaussian beam 3.24 X 1.37 mas in PA —8.9° and epoch C on the right with 
peak flux density 7.7 mjy beam -1 and convolved with a Gaussian beam 2.76 X 1.56 mas in PA —20.5°. 



Epoch 


Orbital 


Total flux density 


Unresolved core 


Jet 


Off source r.m.s. 


(1998) 


phase 


(mjy) 


(mjy beam -1 ) 


(mJy) 


(mjy beam -1 ) 


A 


0.96 


13.0 


6.9 


6.1 


0.04 


B 


0.32 


11.1 


8.0 


3.1 


0.04 


C 


0.68 


11.9 


6.0 


5.9 


0.05 



Table 1. Image properties for Cygnus X-l at 8.4 GHz in August, 1998. Flux calibration is typically accurate to 5 per cent for VLBA 
images at this frequency. The unresolved flux densities per beam were derived from the uniformly weighted (robust -5) images, which 
are not presented in this paper. 



discrete plasmons; i is the angle to the line-of-sight. If the 
source is optically thick in places then the full radiation 
transfer equation must be solved, and then Equation 111 gives 
only an approximation to f3 (Lind & Blandford 1985), see 
e.g. Cawthorne (1991) for a derivation of the (rather small) 
correction factors for a partially optically thick and contin- 
uous jet. 

As we cannot measure the motion of individual compo- 
nents in both the receding and approaching jets we cannot 
directly derive both /3 and i and hence the index k — a. 
However we note that the average accretion disk axis, i, 
(and hence the jet) is likely to be around 40° to the line-of- 
sight (Bruevich et al. 1978) although this is not well defined 
(Brocksopp et al. 1999b). We use the angle to the line-of- 
sight to derive /3 for various values of k — a. 

We assume that we see only one side of the jet in our 
epoch A and use a ratio of approaching and receding flux 
densities of 50 to derive lower limits to the jet speed. The 
ratio is derived by assuming the receding jet flux is < 6 
beams X r.m.s noise per beam, with 6 beams the extent of 



the approaching jet. Although the unresolved flux monitor- 
ing suggests a spectral index of around 0, using -0.6 gives a 
lower limit to (it is unlikely that a synchrotron spectrum 
is much steeper than -0.6). We derive j3 for spectral in- 
dices of and -0.6 and for both the continuous and discrete 
cases (Table 0) . Errors from the poorly constrained angle to 
the line-of-sight dominate these calculations and the results 
should be regarded as order-of-magnitude only, especially as 
any possible attenuation of radiation from the receding jet 
by external free-free absorption cannot be quantified with a 
single frequency. The derived lower limits to J3 suggest the 
jet is at least moderately relativistic > 0.6 c. 

The longer timescale precessional period (142 days) of 
Cygnus X-l, where the deviation of the plane of the disk 
from the orbital plane, 8, may reach as much as ~37° (Brock- 
sopp et al. 1999b) is a further uncertainty. Soft X-ray moni- 
toring (Brocksopp et al. 1999b) suggests that the soft X-ray 
intensity is mid-way between extrema during our observa- 
tions; therefore the disk is close to its average inclination. 

The emitted synchrotron intensity of the jet scales with 
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Figure 3. A high resolution (robust 0) image of Cygnus X-l at 
8.4 GHz; lowest contour 0.15 mjy, convolved with a Gaussian 
beam 2.25 X 0.86 mas in PA -12.4°. 



Spectral index 


Jet type 


lower limit to f3 





continuous 


0.97 





discrete 


0.74 


-0.6 


continuous 


0.83 


-0.6 


discrete 


0.63 



Table 2. Lower limits to the jet velocity as a fraction of c for 
various jet compositions. 



the Doppler factor, D (e.g. Cawthorne 1991), to a power > 
(2 - a). For 13 of 0.75 and i = 40° we find D = 1.55. Using 
the observed long term variations in total flux of < 3 mjy 
gives a maximum value for 8 (assuming, say, a — for the 
radio jet) of only 4°. Therefore either the jet velocity is non- 
relativistic (and hence D ~ 1) or the actual precession of the 
jet is small compared to that suggested for the soft X-ray 
emitting region of the disk. Given our lower limits to /3 (in 
Table 0) we suggest that the inner part of the disk may be 
precessing at significantly less than the angle of around ~37° 
suggested by soft X-ray monitoring. One solution to this 
would be a twisted accretion disk as postulated for SS433 
(Sharp et al. 1984). Alternatively the ~142 day 'period' is 



a timescale for changes in the accretion rate, as is likely in 
GX 339-4 (Corbel et al. 2000). The extended emission on 
VLA scales (Marti et al. 1996) may be consistent with our 
measured jet position angle, particularly if the several mas 
per year south-west proper motion is considered. 

3.3 A continuous and bending jet? 

The higher resolution (robust or partial uniform weight- 
ing) image in Figure K| was made using epoch A, the best set 
of 8.4 GHz data. Even with a resolving beam significantly 
smaller than that provided by natural weighting the jet still 
appears continuous in Figure H, although with small varia- 
tions along it and a general decrease in brightness along the 
jet. Even at this resolution we may find that we are not re- 
solving multiple discrete components (as was not discovered 
in SS433 until detailed VLBI observations were made). 

Epoch A shows a pronounced bend or kink at around 
7 mas from the radio core. The jet region at 8-15 mas sub- 
tends a similar position angle from the core (21-24°) in all 
3 epochs, whereas the inner jet (< 7 mas) is at 17° during 
epoch A. The apparent change in direction of the jet sug- 
gests that the position angle of ejection is varying, perhaps 
either jet precession or jitter, both of which are observed in 
SS433 (Margon & Anderson 1989). It is also possible that the 
kinks observed in the extended emission are consistent with 
a bending jet flow, rather than a ballistic motion, although 
apparently unrelated to the direction of the proper motion 
of Cygnus X-l as would be expected from ram pressure con- 
siderations. Projection effects will enhance the appearance 
of a jet bend (as with apparent 90° bends in blazar jets). 
The higher resolution image in Figure shows that the bend 
occurs in a minimum of surface brightness. This is consis- 
tent with an underlying helical magnetic field in which the 
magnetic field pitch angle reverses sign at the bend in ra- 
dio structure e.g. Laing (1981); Papageorgiou (priv. comm.). 
Measurements of any linear polarization could confirm this 
model. 

The bend is no longer apparent in epochs B and C, 
suggesting that the feature disappears on a timescale of < 
2 days. This could be due to rapid motion and decay via 
adiabatic expansion, or by synchrotron losses. In the latter 
case the magnetic field in the component would need to be 
greater than 16 G, rather higher than is commonly deduced 
for microquasar jets. 

It is also apparent in Figure H that the jet is not resolved 
perpendicular to the flow. Using a beam size of 1 mas across 
the jet and a maximum distance from the core of 15 mas on 
the plane of the sky, the opening semi-angle of the radio- 
emitting material is constrained to <2°. 

The flat radio spectrum in Cygnus X-l specifically 
has been modelled as slowed expansion in a continuous jet 
(Hjellming & Johnston 1988). The predicted spectrum for 
this geometry is flat until a turnover at around 10 GHz for 
an inclination of 40° . However the self-absorption in a slowed 
expanding jet model cannot explain the observed flat spec- 
trum to mm wavelengths (Fender et al. 2000). In comparison 
jet components from GRS 1915+105 are observed to become 
optically thin less than a few minutes after ejection at mm 
wavelengths (Mirabel et al. 1998) . Possibly another inverted 
spectrum component (such as the disk or stellar wind) as 
well as the partially-self-absorbed conical jet is required to 
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explain the wavelength dependence up to mm wavelengths 
of the integrated spectrum. 



4 CONCLUSIONS 

Cygnus X-l has been detected with the VLBA at 15.4 GHz, 
giving a 13 year timebase for proper motion measurements. 
Subsequent imaging at 8.4 GHz shows that a relativistic jet 
with either a bending flow or a variable ejection angle exists 
in Cygnus X-l. The jet appears to be narrow (< 2° opening 
angle). The source was known to be in a low/hard X-ray 
state at the times of all our observations. Given the doc- 
umented X-ray/radio correlations, and our three separate 
images of extended radio emission, it is likely that a jet al- 
ways exists in the 'quiescent' low/hard X-ray state thought 
to harbour a central ADAF. It is unclear how the low-level 
radio flaring of the low/hard X-ray state fits into the picture 
of a steady-state accretion disk input/output system. 

Further VLBI observations with a global VLBI array 
were undertaken in May 2001. A sequence of images within 
a day should show if the radio jet is continuous or consists 
of moving discrete components. 
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